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Abstract. 

In this article we review the standard versions of the Central and of the Levy-Gnedenko Limit 
Theorems, and illustrate their application to the convolution of independent random variables asso- 

— 

ciated with the distribution (X) = 1 + {q - l)3§q{X - figf {jz/g > Q; SSg > Q\ q < 3), 

known as ^-Gaussian. This distribution emerges upon extremisation of the nonadditive entropy 
Sq = k{\ — j[p {xyf dX) I {\~q), basis of nonextensive statistical mechanics. It has a finite vari- 
ance for q < 5/3>, and an infinite one for q> 5/3. We exhibit that, in the case of (standard) inde- 
pendence, the g'-Gaussian has either the Gaussian (if ^ < |) or the a-stable Levy distributions (if 
^ > |) as its attractor in probability space. Moreover, we review a generalisation of the product, the 
<7-product, which plays a central role in the approach of the specially correlated variables emerging 
within the nonextensive theory. 
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INTRODUCTION 

Science, whether in its pure or applied form, is frequently related to the description of 
the behaviour exhibited by systems when some quantity approaches a particular value. 
We might consider the effect on the state of a particle as we furnish it a certain energy, 
E, when this quantity tends to some critical value, Ec, or the response of a system when 
the number of elements goes to infinity, like it usually occurs in thermodynamics and 
statistical mechanics, i.e., the thermodynamic limit. In the latter case, we may focus 
on the outcome of the addition, or arithmetic average, of a large (infinite) sequence 
of random variables associated with a certain observable. This constitutes the basis 
of the celebrated central limit theorem (CLT), which is at the core of the theory of 
probabilities and mathematical statistics. The CLT has its origin at the weak law of 
large numbers of JACOB Bernoulli [1]. For independent random variables, it had its 
first version introduced by Abraham de Moivre in 1733 [2], who used the normal 
distribution to approximate the functional form of a binomial distribution for a large 
number of events. In 1812, his result was later extended by Pierre-Simon Laplace, 
who formulated the now called Theorem ofde Moivre-Laplace [3]. Laplace also used the 
normal distribution in the analysis of errors in experiments, but it was Carl Friedrich 
Gauss who first proved, in 1809, the connection between error in measurement and the 



normal distribution. It is due to this relation that the normal distribution is widely called 
in Physics as Gaussian distribution. Although discovered more than once by different 
people, the fact is that only in 1901, mathematician Aleksandr Lyapunov defined 
the CLT in general terms and proved it in a precisely mathematical fashion [4, 5]. 

After the establishment of a central limit theorem for the addition of independent ran- 
dom variables with finite second-order moment, other versions have appeared, namely, 
the Levy-Gnedenko extension for the sum of independent random variables with diverg- 
ing second-order moment [6, 7], the m-dependent central limit theorem, martingale cen- 
tral limit theorem [8], the central limit theorem for mixing processes among others [9- 
11, 14-20]. 

In this article, we review the fundamental properties of nonadditive entropy, Sq, its 
optimising distribution (known as ^-Gaussian), and the ^-product, a generalisation of the 
product [21, 22] formulated within nonextensive statistical mechanics. We analyse, both 
analytically and numerically, the sum of conventional independent random variables and 
show that, in this case, the attractor in probability space is the Gaussian distribution if 
random variables have a finite second-order moment, or the a-stable Levy distributions 
otherwise. 



NONADDITIVE ENTROPY Sq 

Statistical mechanics, i.e., the application of statistics to large populations whose state 
is governed by some Hamiltonian functional, is strongly attached to the concept of 
entropy originally introduced by RUDOLF JULIUS EMMANUEL Clausius in 1865 [23]. 
The relation between entropy and the number of allowed microscopic states was firstly 
established by LuDWiG Eduard Boltzmann in 1877 when he was studying the 
approach to equilibrium of an ideal gas [24]. Mathematically, this relation is, 

S^kXnW, (1) 

where A: is a positive constant and W the number of microstates compatible with the 
macroscopic state. This equation is known as Boltzmann principle. 

When a system is not isolated, but instead in contact with some kind of reservoir, 

it is possible to derive, from Eq. (1) under some assumptions, the Boltzmann-Gibbs 

w 

entropy, Sbg = ~^ E Pi ^^Pi, where pt is the probability of microscopic configuration 

/ [25]. Boltzmann-Gibbs statistical mechanics is based on the molecular chaos [24] and 
ergodic [26] hypotheses [27]. It has been very successful in the treatment of systems in 
which short spatio/temporal interactions dominate. In this case, ergodicity and indepen- 
dence are justified and Khinchin's approach to Sbg is valid [26]. Therefore, it appears 
as entirely plausible that physical entropies other than the Boltzmann-Gibbs one, can be 
defined in order to treat anomalous systems, for which ergodicity and/or independence 
are not verified. 

Inspired by this kind of systems it was proposed in 1988 [28] the entropy Sq = 
^{ 1~ Lpfj/ll"?) lim^-^i Sq = Sbg) as the basis of a possible extension 



of Boltzmann-Gibbs statistical mechanics [29, 30] where the entropic index q should 
be determined a priori from microscopic dynamics. Just like Sgc^ Sq is nonnegative, 
concave (\/q > 0), experimentally robust (or Lesche-stable [31]) (iq > 0), composable, 
and leads to a finite entropy production per unit time [32, 33]. Moreover, it has been 
shown that it is also extensive [15, 34, 35, 38], hence in compliance with Clausius 
concept on macroscopic entropy and thermodynamics, for a special class of correlated 
systems. More precisely, systems whose phase-space is occupied in a (asymptotically) 
scale-invariant manner. It is upon this kind of correlations that the ^-generalised Central 
Limit Theorems are constructed. 

At this stage let us emphasize the difference between the additivity and extensivity 
concepts for entropy ^ . An entropy is said to be additive if [36] if for two probabilistically 
independent systems, let us say A and B, the total entropy equals the sum of the 
entropies for the two independent systems, i.e., S{A + B) = S{A) +S{B). According 
to this definition, Boltzmann-Gibbs entropy, Sbg^ and Renyi entropy, 5^ [37], 5^ = 
j-r^ln[£"^jp"], are additive, while Sq (q 1), among others [11], is nonadditive. 

Despite the fact of being nonadditive, Sq, just as additive entropies Sbg and S^, is 
composable, as already mentioned. By this we mean that, for a system composed by 
two independent subsystems, A and B, if we know the entropy of each sub-system, then 
we are able to evaluate the entropy of the entire system. Composability for Sbg and S% 
is a consequence of its additivity, whereas for Sq it results from the fact that, considering 
independent subsystems, the total entropy satisfies 

Sq{A+B) _Sq{A) ^ Sq{B) ^ ^^ Sq{A) Sq{B) ^ 

On the other hand, an entropy is defined as extensive whenever the condition, 

iim:^ = .e(0,oo), (3) 

is verified (A'^ represents the number of elements of the system). In this definition, the 
correlation between elements of the system is arbitrary, i.e., it is not important to state 
whether they are independent or not. If the elements of the system are independent (e.g., 
the ideal gas, the ideal paramagnet), then the additive entropies = Sf = Sbg and S^ 
(Va) are extensive, whereas the nonadditive entropy Sq(q^ 1) is nonextensive. In such 
case we have 

5 = 5(1), (4) 

where 5(1) is the entropy of one element when it is considered as isolated. Further- 
more, for short-range interacting systems, i.e., whose elements are only asymptotically 
independent (e.g., air molecules at normal conditions), i.e., strict independence is now 

violated, Sbg and 5^ still are extensive, whereas Sq (q 7^) still is nonextensive. For such 
systems, sbg = limA?^oo5BG(^)/^ 7^ "^(1)- Conversely, for subsystems of systems ex- 
hibiting long-range correlations [12], Sbg and 5^ are nonextensive, whereas Sq can be 



In this discussion we will treat elements of a system as strictly identical and distinguishable. 



extensive for an appropriate value of the entropic index q^l. This class of systems has 
been coined as q-describable [13]. 



Optimising Sq 

Let us consider the continuous version of the nonadditive entropy Sq, i.e., 

\-f\p{X)]^dX 



Sq — k- 



l-q 



(5) 



The natural constraints in the maximisation of (5) are (hereinafter k= I), J p {X) dX ■ 
1 , corresponding to normalisation, and 



x^Pi^9^dx-(x) -a 



J[pix)]Ux 

corresponding to the q-generalised mean and variance of X, respectively. 
From the variational problem we obtain 



(6) 
(7) 



1-9 
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%{X)^^q[l + {q-l)3§q{X-jlqy 

(if the quantity within brackets is nonnegative, and zero otherwise) where. 
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and = [(3 — ^) o^~\ \ Standard and generalised variances, Oq and <7^ are related 



through = (y^^2^, for ^ < § . 

Defining the q-exponential function ^ as 

e^q^[l + {\-q)x]^^ {e\^e^), 
(e^ = if 1 + ( 1 — q)x < 0) we can rewrite PDF (8) as 



\x) = a/qeq 



(10) 



(11) 



^ Other generalisations for the exponential function can be found at Ref. [39]. 



hereon referred to as q-Gaussian. The inverse function of the ^-exponential, the q- 
logarithm, is In^ {x) = ^^4^ (x > 0) . 

For q = the (/-Gaussian recovers Student's ^-distribution with m degrees of 
freedom (m = 1,2,3,...) and finite moment up to order m. So, for ^ > 1, PDF (11) 
presents an asymptotic power-law behaviour. Complementarily, if ^ = with n = 
3,4,5, p (x) recovers the r-distribution with n degrees of freedom. Consistently, for 

q< I, p{x) has a compact support defined by the condition |x — jU^ | < ^jz^o^ . 



^-calculus 

The nonadditivity property of Sq, assuming/or independent systems A and B, 

Sq {A+B) _ Sq (A) ^ 5, {B) ^ Sq (A) 5, (B) ^ ^^^^ 

has inspired the introduction of a new algebra [21, 22] composed by q-sum, x®qy = 
x + y + {I — q)xy, and the q-product 

x^qy = [x^-* - 1] ^ . (13) 

The corresponding inverse operations are the q-difference,xQqy, and the q-division, x 0q 
y, such that, {x 0qy) 0qy = x. The ^-product can be written by using the basic function 
of nonextensive formalism, the (/-exponential, and its inverse, the ^-logarithm. Hence, 
x0qy = exp^ [in^ x + In^ , which for ^ — > 1 , recovers the usual property In (xxy) = 
\nx + \ny (x^y > 0), where x xy = x®iy. Since exp^ [x] is a non-negative function, the 
^-product must be restricted to values of x and y that respect condition 

\x\^~^ + \y\^~^-l>0 (14) 
We can enlarge the domain of the ^-product to negative values of x and y by writing it as 

x0qy= sign (xy) exp^ [in^ |x | -h In^ | ] . (15) 

We list now a set of properties of the ^-product: 

1. x<S)iy = xy ; 

2. x<Siqy = y0qX ; 

3. {x^qy)0qz^x0q{y0qz) ^x®qy®qZ^ [x^-^ + j^"* + z^-^ - 3] ^ ; 

4. {x®q 1) = x ; 

5. lnq[x®qy] = \nqx + \nqy\ 

6. In^ (xy) = In, (x) + In, (y) + (1 - 1% (x) In, {y) ; 

7. (x(g),j)"^ =x~^ (8)2-9 



8. {x®qO) = 




if (^> 1 andjc>0)orif < 1 and0<jc< 1) , 
if q < 1 andx > 1. 



For special values of q, e.g., q = 1/2, the argument of the (^-product can attain 

nonpositive values, specifically at points for which \x\ ^ + |j| ~^ — 1 < 0. In these 
cases, and consistently with the cut-off for the ^-exponential we have set x®qy = Q. 
With regard to the ^-product domain, and restricting our analysis of Eq. (14) to x,y > 0, 
we observe that for q — oo the region {0<jc< 1,0 <y < 1} leads to a vanishing q- 
product. As the value of q increases, the forbidden region decreases its area, and when 
^ = we have the limiting line given by x+y = 1, for which x(8)oy = 0. Only for ^ = 1, 
the whole set of real values of x and y has a defined value for the ^-product. For q>\, 
condition (14) yields a curve, |x|^^ + |y|^^ = 1, at which the ^-product diverges. This 
undefined region increases as q goes to infinity. At ihcq^oo limit, the ^-product is only 
defined in {x> l,y< l}U{0<x< \,0<y < l}U{jc< l,y> 1}. This entire scenario 
is depicted on the panels of Fig. 1 . The profiles presented by x ®ooy and x®-ooy illustrate 
the above property (8). To illustrate the ^-product in another simple form, we show, in 
Fig. 2, a representation oix®qX for typical values of q. 



LEVY DISTRIBUTIONS 



In the context of the CLT for independent variables, apart from the Gaussian distribution, 

^ iX), another stable distribution plays a key role, the a-Levy distribution, {X) 
(0 < a < 2). If 5f {X) is characterised by its fast decay, {X) is characterised by its 
'fat' tails, since it allows both small and large values of X to be effectively measurable. 
Widely applied in several areas, (X) is defined through its Fourier Transform [6], 

Ll{k) = exp -a|)t|"|l + /fitan(af)^} , (ay^l), (16) 

(where a is the Levy exponent, and B represents the asymmetry parameter). By this we 
mean that Levy distributions have no analytical form in X, excepting for special values 

of a = 1. For B = 0, L^^= La, the distribution is symmetric. Regarding a values, one 
can verify that Li ix) corresponds to the Levy-Smirnov distribution, and that Li (X) is 

the Cauchy or Lorentz distribution (Li (X) coincides with the ^2 (X) distribution). For 
a = 2, La (k) has a Gaussian form, thus the corresponding distribution is a Gaussian. 
Carrying out the inverse Fourier Transform on La {k),La (X) = J-^ La (k) dk, 

we can straightforwardly evaluate the limit for small X, La {X) ^ ^;raa"/°^^ F [^] , 
and the limit for large X, 

, , aar[a] sinful 

^^TT^, (17) 

From condition, < a < 2, it is easy to prove that variables associated with a Levy 
distribution do not have a finite second-order moment, just like '^q iX) with \ <q<'i 



FIGURE 1. Representation of the ^-product, Eq. (15), for q = -oo, -5, -2/3, 0, 1/4, 1, 2, oo. As it is 
visible, the squared region {0<x< l,0<y< l}is gradually integrated into the nontrivial domain as q 
increases up to ^ = 1 . From this value on, a new prohibited region appears, but this time coming from large 
values of ( |x| , IjI ) . This region reaches its maximum when g = oo. In this case, the domain is composed by 
a horizontal and vertical strip of width 1. 



(see [40] and references therein). In spite of the fact that we can write two distributions, 
(X) and La {X), which present the same asymptotic power-law decay (with a = |zf )> 
there are interesting differences between them. The first one is that, as we shall see later 



0.0 0.5 1.0 ^1.5 2.0 



FIGURE 2. Representation of the ^'-product, x(^qX for q = —5, 0, 1,2, °o. Excluding q ~ I, there 
is a special value x* — 2'/^''^'', for which q < I represents the lower bound [in figure x* {q ~ —5) = 
2-1/^ ~ 0.89089 and x* 0) = 1 /2], and for ^ > 1 the upper bound [in figure x* {q^2)^ 2]. For 
q = ±oo, x®qX lies on the diagonal of bisection, but following the lower and upper limits mentioned 
above. 



on, La {X) together with 5f {X) are the only two stable functional forms whenever we 
convolute independent variables [20]. Therefore, distribution ^g(X) {q ^ 1,2) is not 
stable for this case. The other point concerns their representation in a log — log scale. 
Contrarily to what happens in log — log representations of ^q{X), an inflexion point 
exists at Xi for L« {X) if 1 < a < 2, see Fig. 3. Since in many cases the numerical ad- 
justment of several experimental/computational probability density functions for either 
a Levy or a ^-Gaussian distribution seems to be plausible, the presence of an inflexion 
point might be used as an extra criterion to conclude which one is the most adequate. 
This has clear implications on phenomena modelling. 




FIGURE 3. Left panel: Gaussian and a-stable Levy distributions for a approaching 2 in Eq. (16) with 
a = \ and B = 0. As referred in the text for values of a closer to 2, Levy distribution becomes almost 
equal to a Gaussian up to some critical value for which the power law behaviour emerges. Right panel: 
Locus of the inflexion point of the a-stable Levy distributions, Eq. (16), with a = 1 and B = Q. Conti-arily 
to what happens with {X), when Levy distributions are represented in a log-log scale, they exhibit an 
inflexion point which goes to infinity as a ^ 1 (Cauchy-Lorentz distribution (^)) and a 2 (Gaussian 
distribution ^ (Z)) too. We also show the projections onto the planes —Xj, — a, and a—Xi. 



CENTRAL LIMIT THEOREMS FOR INDEPENDENT VARIABLES 



Variables with finite variance 



Let us consider a sequence Xi , X2, . . . ^X^ of random variables which are defined on the 
same probability space, share the same probability density function, p (X), with cr < 00, 
and are independent in the sense that the joint probability density function of any two Xi 
and Xj, P {Xi, Xj) , is just p (Xi) p (Xj) . 

Hence [41, 42], a new variable 



X1+X2 + ...+XN 



N 



(18) 



A' 



with raw moments, {Y") = (^-^ j y > has its probability density function given 

by the convolution of N probability density functions, or, since variables are indepen- 
dent, by the Fourier Transform, 



-^[PiYMk) = l{le^'"^p{X)dx] 

= ^ exp \n In [1 + ik^^ - i^2i^ + O (Ar-3) 



N 



(19) 



Since ln(l+x)=jc— y + y + ..., expanding up to second order in X we asymptotically 
have, in the limit N ^ 0° , 



^[P{Y)]{k) ^ ^exp 



NUik^ 



N 2N^ 



[X') - {xy 



(20) 



Defining jix = {X) , as the mean value, and = (X^) — (X)^ , as the standard deviation 
we have 

'^[P{Y)]{k) ^ ^cxp[-ik^x-if(y^] . (21) 

Performing the Inverse Fourier Transform for J^[P{Y)] (k), we obtain the distribution 
PiY), 

p{Y) = Je-'''^^[P{Y)]{k)dk , (22) 

which yields, 

lg-(i'-Mz)V(2#) 



P{Y) 



(23) 



Remembering that, from Eq. (18) jUy = /i^ and Oy — OxN we finally get 

P(y) = _L_g-(i'-|ty)V(2<T2) ^ (24) 

which is a Gaussian. It is also easy to verify that P (^Z — — \/NP{Y), i.e., the 



Gaussian distribution scales as A^^ 2 . 



Example: The convolution of independent q-variables 



Let us assume that the probabiUty density function (X) is a (^-Gaussian distribution, 
^q{X), which, as a simple illustration, has q=\, and <7 = 1 < oo ^. In this case, the 

Fourier Transform for ^3 (X) is ^ ^3 (X) {k) = (1 + exp [- \k\\. The distribution, 

^ (F), where Y = Xx +X2 + . . . +Xiv, is given by 



^{Y) = ^Jcxp[-ikY][^[^3{X)\{k)Ydk . 



(25) 



1 



Expanding |^ [^3 (X) (/:)|^ around A: = we obtain, |^ (X) (/:)|^ ~ 
^Nk^ + hN\k\^. From the CLT, we observe that distribution ^ (F), for large N, is well 



2i»n, -r 3 

described by a Gaussian, ^ (F) 



exp 



2Af 



, at its central region (for large). 



Because of singularity \k\ , associated with the divergence in n-order statistical moments 
(n > 3), we have the remaining distribution described by a power-law, ^ (F) ~ ^Y^, 
for large F. This specific behaviour is depicted in Fig. 4. Therein we observe a crossover 
from Gaussian to power law at YN~^/^ of order \/\nN, which tends to infinity. 
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FIGURE 4. Both panels represent probability density function ^{Y) vs. Y (properly scaled) in log- 
linear (left) and log-log (right) scales, where Y represents the sum of independent variables X having 
a {X) distribution. Since variables are independent and its variance is finite, £P (Y) converges to a 
Gaussian as it is visible. It is also visible in the log-linear representation that, although the central part of 
the distribution approaches a Gaussian, the power-law decay subsists even for large as it is depicted in 
log-log representation. 



^ The value ^ = 5 appears in a wide range of phenomena which goes from long-range hamiltonian systems 
[43^5] to economical systems [46]. 



Variables with infinite variance 



Figure 4 is similar to what is presented in Fig. 3 (left panel). As it is visible there, as 
a goes to 2, the distribution La {X) nearly collapses onto the Gaussian distribution up to 

some critical value X* . Beyond that point, the asymptotic power-law character emerges 

™ 1 

and the distribution falls as |X | 

If we consider the sum, F = Zi +Z2 + . . . +^7^, of A^^ random variables, X\ , X2, . . . , Xn, 
which share the same Levy distribution, La{X). The distribution ^(Y)is then given by, 

^(Y) = i^I^:e-"' {La{k)fdk 

= ^J+::cxp[-ikY-aN\k\'']dk. ^ ^ 

Introducing a new variable, co = kN^/", we get ^ (Y) = N~^/"La (^}^^ j • Hence, for 

Levy distributions we also have the scaling property, but with exponent I /a. This gen- 
eralised version of the Central Limit Theorem, originally due to Gnedenko-Kolmogorov 
[7], is applied to any distribution p{x) which, in the limit x — * 0°, behaves as 
(/i < 2). In other words, the probability density function of the sum of N variables, each 
one with the same distribution p{x) ~ (0 < }! < 2), converges, in the N ^ 00 

limit, to a a-stable Levy distribution. 



FINAL REMARKS 

In this article we have reviewed central limit theorems for the sum of independent 
random variables. As we have illustrated, in the absence of correlations, the convolution 
of identical probability density functions which maximise non-additive entropy Sq, 
behave in the same way as any other distribution. In other words, when the entropic 
index ^ < | the variance of {X) is finite, hence the convolution leads to a Gaussian 
distribution. On the other hand, i.e., when ^7 > | the variance diverges. As a consequence 
the convolution leads to a a-stable Levy distribution with the same asymptotic power- 
law decay of '^q {X) (the marginal case q = 5/3 yields a Gaussian distribution, but with a 
logarithmic correction on the standard oc t scaling, i.e., it yields anomalous diffusion). 
We have also exhibited that there is an important difference between ^-Gaussian and 
oc-stable Levy distributions, namely the emergence of an inflexion point on the latter 
type when they are represented in a log-log scale. In the subsequent paper (Part 11) we 
will show that the strong violation of the independence condition introduces a drastic 
change in the probability space attractor. Specifically, for a special class of correlations 
(q-independence), it is the ^-Gaussian distribution which is a stable one. 



ACKNOWLEDGMENTS 

We are deeply thankful to E.P. Borges and R. Hilfer for fruitful discussions related to 
both Part I and Part II. Partial financial support from Pronex, CNPq, Faperj (Brazilian 
agencies) and FCT/MCES (Portuguese agency) is acknowledged as well. 



REFERENCES 



1. J. Bernoulli, Ars Conjectandi (Basel, 1713). 

2. A. de Moivre, The Doctrine of Chances (Chelsea, New York, 1967). 

3. P.-S. Laplace, Theorie Analytique des Probabilites (Dover, New York 1952). 

4. A. M. Lyapunov, Bull. Acad. Sci. St. Petersburg 12 (8), 1 (1901). 

5. H. Tijms, Understanding Probability: Chance Rules in Everyday Life (Cambridge University Press, 
Cambridge, 2004). 

6. P. Levy, Theorie de I'addition des variables aleatoires (Gauthierr-Villards, Paris, 1954). 

7. B. V. Gnedenko, Usp. Mat. Nauk 10, 1 15 (1944) (Translation nr. 45, Am. Math. Soc, Providence). 

8. P. Hall and C. C. Heyde, Martingale Limit Theory and Its Application (Academic Press, New York, 
1980). 

9. S. Artstein, K. Ball, F. Barthe and A. Naor, J. Amer. Math. Soc. 17, 975 (2004). 

10. C. TsalHs, Quimica Nova 17, 468 (1994). 

11. M. Gell-Mann and C. Tsallis, eds., Nonextensive Entropy - Interdisciplinary Applications (Oxford 
University Press, New York, 2004). 

12. F. Caruso and C. TsalUs, in this same volume. 

13. C. Tsallis, M. Gell-Mann and Y. Sato, in Nonextensive Statistical Mechanics: New Trends, New 
Perspectives, J.P.Boon and C. Tsallis (eds.), Europhysics News 36 (6) (2005). 

14. C. Beck and F. Schlogel, Thermodynamics of Chaotic Systems: An Introduction (Cambridge Univer- 
sity Press, Cambridge, 1993). 

15. C. TsalUs, M. Gell-Mann and Y. Sato, Proc. Natl. Acad. Sc. USA 102, 15377 (2005); idem, Euro- 
physics News 36, 186 (2005). 

16. L.G. Moyano, C. TsalHs and M. Gell-Mann, Europhys. Lett. 73, 813 (2006). 

17. S. Umarov, C. Tsallis and S. Steinberg, arXiv : cond-mat / 0603593 (pre-print, 2006). 

18. S. Umarov, C. Tsallis, M. Gell-Mann and S. Steinberg, arXiv: cond-mat/06006038 (pre-print, 
2006). 

19. S. Umarov, C. TsalUs, M. Gell-Mann and S. Steinberg, arXiv : cond-mat / 06006040 (pre-print, 
2006). 

20. A. Araujo and E. Guine, The Central Limit Theorem for Real and Banach Valued Random Variables 
(John Wiley & Sons, New York, 1980). 

21. L. Nivanen, A. Le Mehaute and Q.A. Wang, Rep. Math. Phys. 52, 437 (2003). 

22. E.P Borges, Physica A 340, 95 (2004). 

23. E. Fermi, Thermodynamics, (Doubleday, New York, 1936). 

24. L. Boltzmann, Lectures on Gas Theory, (Dover, New York, 1995). 

25. K. Huang, Statistical Mechanics, (John Wiley & Sons, New York, 1963). 

26. A.I. Khinchin, Mathematical Foundations of Information Theory (Dover, New York, 1957) and 
Mathematical Foundations of Satistical Mechanics (Dover, New York, 1960). 

27. E.G.D. Cohen, Pramana - J. Phys. 64, 635 (2005) [Boltzmann Award lecture]. 

28. C. Tsallis, J. Stat. Phys. 52, 479 (1988). 

29. E.M.F Curado and C. TsalUs, J. Phys. A 24, L69 (1991); Corrigenda: 24, 3187 (1991) and 25, 1019 

(1992). 

30. C. TsaUis , R.S. Mendes and A.R. Plastino, Physica A 261, 534 (1998). 

31. B. Lesche, J. Stat. Phys. 27, 419 (1982). 

32. V. Latora and M. Baranger, Phys. Rev. Lett. 273, 97 (1999). 

33. C. Tsallis in Nonextensive Entropy - Interdisciplinary Applications, M. Gell-Mann and C. Tsallis 
(eds.) (Oxford University Press, New York, 2004). 

34. C. Tsallis, Proceedings of the 31st Workshop of the International School of Solid State Physics 
"Complexity, Metastability and Nonextensivity", held at the Ettore Majorana Foundation and Centre 
for Scientific Culture (Erice, July 2004), eds. C. Beck, A. Rapisarda and C. TsaUis (World Scientific, 
Singapore, 2005); Y. Sato and C. Tsallis, Proceedings of the Summer School and Conference on 
Complexity in Science and Society (Patras and Ancient Olympia, 14-26 July 2004), Complexity: An 
unifying direction in science, eds. T. Bountis, G. Casati and I. Procaccia , Int. J. Bifurcation and Chaos 
16, 1727 (2006). 

35. J. Marsh and S. Earl, Phys. Lett. A 349, 146 (2005). 

36. O. Penrose, Foundations of Statistical Mechanics, (Pergamon Press, Oxford, 1970). 



37. A. Renyi, Proceedings of the 4th Berkeley Symposium on Mathematics, Statistics and Probability, 

547(1960). 

38. J. A. Marsh, M.A. Fuentes, L.G. Moyano and C. Tsallis, Physica A 372, 183 (2006). 

39. T. Ernst, J. Nonlinear Math. Phys. 10, 487525 (2003). 

40. D. Prato, C. Tsallis, Phys. Rev. E 60, 2398 (1999). 

41. J.W. Lindeberg, Math. Z. 15, 211 (1922). 

42. W. Feller, An Introduction To Probability Theory And Its Applications, V. I (John Wiley & Sons, New 
York, 1968). 

43. V. Latora, A. Rapisarda and C. TsalUs, Phys. Rev. E 64, 056134 (2001). 

44. C. Tsallis, Physica A 344, 718 (2004). 

45. L.G. Moyano, C. Anteneodo, Phys. Rev. E 74, 021 1 18 (2006). 

46. R. Osorio, L. Borland and C. Tsallis, in Nonextensive Entropy - Interdisciplinary Applications, M. 
Gell-Mann and C. Tsallis (eds.) (Oxford University Press, New York, 2004); S. M. Duarte Queiros, 
Quantitat. Finance 5, 475 (2005). 



